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Abstract 

The statistics of quantum Poincare recurrences in Hilbert space for dia- 
magnetic hydrogen atom in strong magnetic field has been investigated. It 
has been shown that quantities characterizing classical chaos are in a good 
agreement with the ones that are used to describe quantum chaos. The 
equality of classical and quantum Poincare recurrences has been shown. 
It has been proved that one of the signs of the emergence of quantum 
chaos is the irreversible transition from a pure quantum mechanical state 
to the mixed one. 

PACS:05.45.Mt 



1 Introduction 

The investigation of the chaotic processes in the deterministic systems is one of 
the most important trends of the modern physics [1-5]. The successes of this 
current for nonlinear classical systems that are described via trajectories of the 
system in Hilbert space are obvious [6-8] . Yet the conception of deterministic 
chaos in quantum systems, where the idea of trajectory is inapplicable, stays a 
subject for serious discussions so far. There is no common understanding of the 
problem of quantum chaos for today. Usually the quantum chaos is considered 
as a set of phenomena taking place in quantum systems, the classical analogues 
of which display chaos [9-16]. 

Generally when investigating the problems of quantum chaos, one has to 
deal with the Hamiltonians of the form: 

H{p,q,\)^Ho{p,q) + \V{q), (1.1) 

where p, q are a set of classical coordinates and momentum, Hq {p, q) is classically 
integrable part of Hamiltonian, V(q) is a part of Hamiltonian which leads to 
the nonintegrability of the classical equations of motion when added, A is a 
parameter by variation of which the system may be driven in the domain of 
chaotic dynamics. 
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The goal of this work is the investigation of a concrete physical system 
described by the Hamiltonian of the form (1.1). The object of our investigation 
is a diamagnetic hydrogen atom placed in a strong magnetic field. Diamagnetic 
hydrogen atom, subjected to the influence of the strong magnetic field, has been 
considered in the quasi-classical approximation in the ref [16]. Namely a photo- 
absorbtion spectrum in the transition regime to chaos has been studied. The 
main purpose of our work is the investigation of the irreversible evolution of 
the quantum-mechanical system and the passage from the quantum-mechanical 
to the quantum-statistical consideration. We shall try to compare and find out 
the link between the phenomena taking place under both classical and quantum 
considerations. 

Relatively weak magnetic field applied to atom leads just to the energy levels 
splitting, i.e. to the Zecman effect. In case when the hydrogen atom is placed 
in the extreme high magnetic field, the motion of the electron in the plane 
perpendicular to the field direction is completely defined by this magnetic field 
and not by the coulomb field of the nucleus, and because of this, the atom is 
found to be deformed in the transverse direction. At the same time the field 
does not aff'ect on the lengthwise motion, and in this direction the atom keeps 
its size. 

One may easily estimate the value of the magnetizing force when the defor- 
mation of the electron shell takes place. One has to compare the value of Bohr 
radius ag = fi^/mge^ (where h is a Plank constant, me and e are the electron 
mass and charge respectively) with a characteristic size of the localization do- 
main of electrons in the magnetic field in the ground state when n = 0, s = 0, 
as = ^hc/eB. If < oo, then the magnetic field has the defining effect. This 
condition leads to the following estimate of the extremely strong field 
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B > ^^^^ = Be = 2.35 • lO^Ts, (1.2) 

in which the atom is deformed. Notice, that the magnetic fields required for 
this, the intensity of which is defined by the condition (1.2), according to the 
modern concepts, may exist on the surface of some astrophysical objects. The 
neutron stars (or pulsars) appearing as the result of the collapse of the super- 
nova, belong to them. That is why the investigation of the diamagnetic atoms 
in the conditions of extremely high magnetic fields is of doubtless interest. 

2 Classical Consideration 

Consider the hydrogen atom placed in the magnetic field directed along z-axis 
B = Bcz ■ We shall use Larmor theorem. According to this theorem, for study- 
ing the diamagnetism theory for single atoms it is enough to consider the motion 
of the electron in the coordinate system, rotating with the Larmor frequency 
ojl = "feB, where is a gyromagnetic ratio for electrons. 

During the transition to the rotating coordinate system, the change of the 
kinetic energy of the electron will occur. This change is conditioned by the 
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addition to the velocity of in the rotating coordinate frame u the cross-product 
of the Larmor frequency ul and the vector r of the location of electron. Taking 
into account, that the direction of vector Ql is opposite to the vector 7e-B, and 
gyromagnetic ratio of the electron is 7,; = — ^ (here c is the velocity of light) , 
we shall achieve for the velocity of electron in a motionless coordinate frame; 

6 -* 

V = u-\- [uJLr\ =u + \Br\. 

ZnieC 

Then the kinetic energy of the electron will be of form: 

K = ^meU^ + me[iOL'r]u + ^meWcuLrW^ = 

One may easily see, that if z-projection of the angular momentum is equal 
to zero, L2 = the task can be reduced to the diamagnetic Kepler problem. 
Finally in the system of atomic units the Hamiltonian of the considered problem 
takes the form: 

H=\p^-- + \\{x^+y^) = E, (2.1) 
2 r 8 

where A = " a -B is a constant describing the connection of the system with 

m^ce ^ ^ 

the external magnetic field. The form of the Hamiltonian (2.1) is similar to the 
one of (1.1). For A = the equations corresponding to (2.1) are integrable. For 
A ^ the system is nonintegrable. Finally with the increase of the amplitude 
of the magnetic field, starting with some definite value of Ac, the dynamics 
appears to be chaotic. The use of the numerical methods for the Hamiltonian of 
the form (2.1) is connected with definite diSiculties because of the appearance 
of singularity in the r = point. 

To get over this problem it is necessary to use scmiparabolic coordinates 
that are connected with the Cartesian ones via formulas: 

= \/r + z 
V = y/r - z 

and make a transition to a non-homogeneous time via transformation dt = 
(i/^(t) -|- ii'^{T))dT [15]. Then taking into account that = 0, the Hamiltonian 
(2.1) can be transformed to the following form [15]: 

\{pI+pI) - Eif^' + .2) + iAM^.2(M^ + = 2. (2.2) 

I.e. the Hamiltonian is of form 

H{p^,p„,H,v,X) = Ho{p^i,Pt.,lJ-,iy) + XV{tJ,,u), (2.2a) 

where Hq = ^{pfi +Pu) ~ -^'(a*^ + '^^) is the integrable part of the Hamiltonian; 
V = g/i^z^^(/ti^ + u'^) is the nonintegrable part of the Hamiltonian. 
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Figure 1: Phase trajectory on the piano of variables {p^,fi), obtained after the 
numerical integration of the Hamilton equations, corresponding to the Hamiltonian 
(2.2). The result (as all the other numerical calculations given in this work) is 
achieved by using the software MathWorks MatLab 6.0 for the values of parame- 
ters E = —0.125 ~ 10^^^ J, A — 1.0. As it is seen from the graph, the motion is of a 
chaotic type, and the destruction of the invariant tori points to this. 



As one may see from (2.2a), the transformed Hamiltonian by its form co- 
incides with the Hamiltonian of interacting oscillators, and in (2.2) the singu- 
larities which are typical for the Hamiltonian (2.1) are absent. Though one 
should take into account, that during the direct numerical calculations, right 
after the integration of the Hamilton equations corresponding to (2.2a), one has 
to perform a reverse transition to a real time. 

Henceforth the transformed Hamiltonian will be the object of our investiga- 
tions in case of classical consideration. 

As it is clear from (2.2), the value of the magnetic field B is the parameter of 
nonlinearity. We shall study a case of a strong nonlinearity, i.e. the case when 
the methods of perturbation theory are inapplicable. In such case the chaos 
may appear in the system. We should expect to get a chaos in case of high 
magnetic fields B > Be- 
When chaos appears, the mechanical concepts lose their meaning. For de- 
scribing the processes taking place in the system, the statistical concepts: Kol- 
mogorov entropy [7] and the fractal dimension of the phase space of the sys- 
tem [17-19] become significant. The study of the statistical features of the 
dynamical system can be performed only numerically. 

The results of the numerical computations, gained by solving the Hamilton 
equations, corresponding to (2.2), are given in Fig.l and Fig. 2. 

As seen from Fig.l and Fig. 2, the destruction of the invariant tori took 
place. The motion is of the chaotic type. For more convincing proof of the 
chaotic character of motion one has to perform Fourier analysis of the obtained 
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Figure 2: Phase trajectory on the plane of variables {pi^,v) plotted under the condi- 
tions similar to the one in Fig.l. The destructed invariant torus, according to KAM 
theorem [7] , testifies to the existence of the classical chaos. 



numerical data. 

To define Kolmogorov entropy one has to perform a Fourier transformation 
of tlie correlation function 



G.(T) = Mi + T)Ki)), 
G(w) = / G{T)e"^^dT ^' 



(2.3) 



1 + a;2r2 



where 

((•••))= lim ^ f{---)dt, 

T^oo 1 Jo 

means averaging over the time, Tc is the correlation time, connected with the 
Kolmogorov entropy by the ratio ho ^ ^. In practical calculations of (2.3) it 
is more convenient to use a realization of a Fast Fourier Transform algorithm 
[20,21]. As the numerical calculations show, the value of the correlation time is 
Tc ~ 10^^''^s. The results of numerical calculations are presented in Fig. 3. 

There is one more interesting property of dynamical systems connected with 
the Kolmogorov entropy that describes chaotic motion and is called Poincare 
recurrence. The main point of the Poincare recurrence is that any system with 
a finite phase space after some definite time returns to its initial state [22] . The 
study of the time of Poincare recurrence is one of the most powerful methods 
of the analysis of the nonlinear systems. 

It turned out to be that in conditions of strong chaos (i.e. when the phase 
plane doesn't contain the stability islands) the distribution of the times of 
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Figure 3: The dependence of the Fourier image of the correlation function G^(aj) on 
the frequency w [lO^^Hz), plotted using Eq. (2.3) and FFT method. The presence of 
the finite width of the correlation function Gij,{u>) is the evidence of the appearance of 
chaos. In other words, one may say, that oscillatory process is not characterized by a 
definite frequency (period). All the frequencies (periods) in a finite frequency interval 
5u> ~ 1/tc = 200THz are involved into oscillatory process. 

Poincare recurrences has a Lorenz structure [22] : 

P(t) = — exp(-T/r^ec), (2.4) 

Tree 

where Tree = /o°° TP{T)dT is the average recurrence time connected with Kol- 
mogorov entropy as Tree = Thus, calculating the value of the Kolmogorov 
entropy one can define the distribution of the recurrence times. 

We in details touched upon the question about classical Poincare recurrences, 
because henceforth we will be interested in quantum Poincare recurrences. The 
idea of quantum Poincare recurrence was introduced in the recently published 
work by G. M. Zaslavsky and A. lomin [22]. We shall return to the description 
of this phenomenon during a quantum consideration. 

One of the signs of the appearance of chaos is a discrete (fractal) dimension 
of the phase space of the system. 

For estimation of the fractal dimension we shall use the Grassberger-Procaccia 
algorithm [17-19]. The main point of this algorithm consists into following: 

Assume we have obtained from a numerical solution of the equations of mo- 
tion a set of state vectors {a;,, i = 1,2, - ■ ■ , N} that correspond to the sequential 
steps of integration. In our case Xi denote a complete set of variables, character- 
izing the phase space of the system: Pij,{ti), PuiU), IJ.{ti), p{ti), ti = 0, tff = T, 
where U € [0, T] is the time interval for the numerical integration. Defining 
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Figure 4: Graph of the dependence of C(e) on e. The graph is plotted via the 

integration of classical Hamilton equations corresponding to the Hamiltonian (2.2a) 
for the values of parameters E = —0.125 ^ 10~^^ J, A = 1.0. 

some (small) e, we can use our set for estimation of the following sum: 



According to Grassbcrgcr-Procaccia algorithm, if one knows C(e), the fractal 
dimension of the strange attractor may be defined 



One should calculate C{e) for different values of e and represent the results in 
coordinates log(e) and log(C(e)). 

The expected dependence of C{e) has a form of e^, so that the obtained 
graph must have a straight-line form with D angular coefficient. Results of the 
numerical calculations are presented in Fig. 4 and Fig. 5. 

3 Quantum Consideration 

As it has been already mentioned, usually the quantum chaos is considered as 
a set of phenomena that appear in quantum systems, the classical analogues of 




where 6 is the Heaviside step function 
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Figure 5: Graph of the dependence of C(e) on e, plotted in a logarithmic scale using 
the Grassberger-Procaccia algorithm [17-19]. The graph is plotted via the integration 
of classical Hamilton equations corresponding to the Hamiltonian (2.2) for the values 
of parameters E — —0.125 ~ 10^^^ J, A = 1.0. The curve lino corresponds to the 
results of the numeric data. The straight-line is plotted via the processing of the 
numerical data with least squares algorithm. From this graph one can estimate the 
fractal dimension of the attractor using the formula D = '°g(C(^2))-i°g(C(^i)) ^ 2.857. 

° log(e2)-log(ei) 



which display chaos. The difficulty in definition of quantum chaos is connected 

with a fact, that the concepts characterizing classical chaos (local instability 
of phase trajectories) lose their sense in quantum consideration. Therefore, 
for describing quantum chaos one has to use methods and concepts which are 
typical for quantum mechanics. 

Our goal in case of quantum consideration, is the demonstration of the fact 
that one of the signs of the quantum chaos is a formation of the mixed state. 
For endorsement of this assumption we shall use methods of random matrix 
theory [1,28,29]. Due to this fact, we have to use basis of eigenfunctions of 
the integrable part of the Hamiltonian (2.1) and we cannot expend the wave 
function in Landau states, as was done in [15]. This makes situation more 
complicated but as we are interested in the state I = 0,m = , some analytical 
results still can be obtained. 

Let us write down the Schrodinger equation for the Hamiltonian (2.1) 

H{p,f,X)\i>{f,X)) = E{X)\i,{f,X)). (3.1) 

It is easily seen from (3.1), that the Schrodinger equation in contrast to the 

classical equations of Hamilton is linear. Like in the case of classical systems, 
this conclusion means that Eq.(3.1) is completely deterministic and does not 
contain the properties of chaos, the main sign of which is fortuity and abrupt- 
ness of the behavior of system. Because of this the question appears: Then 
how can one observe nonregular behavior and quantum chaos in the quantum 
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systems? To answer this question we shall use some well-known facts from 
quantum mechanics [23,24,25]. 

Let us define via \'ipn{r,X)) the required eigenfunction of the Eq.(3.1) and 
expand it in eigenfunctions of the Hamiltonian H{p, f, 0) = Ho{p, r) 

IV'n(r,A)) = ^C„„(A)|$„(r)). (3.2) 

m 

Here |$m(r)) is the eigenfunction of the Hamiltonian Ho{p,r). 

According to the universally recognized rules of quantum mechanics, the co- 
efficients of the expansion Cm„(A) taken squared by module, define the proba- 
bility to find the considered system in state described by wave-function \ ^m{r)), 
and that is why they satisfy the conditions of normalization: 

$^|C„„(A)p = l. (3.3) 

m 

Substituting the expansion (3.2) into Schrodinger equation (3.1), one can obtain 
a matrix equation for Cmn 

HkpCpn = EnCknW, (3.4) 

p 

where the elements of matrix Hkp are defined by the integrals of the type: 

Hkp = {<^k{f^\H{p,X)\<^p{f)). (3.5) 

It is shown in quantum mechanics [23-25] that the scries (3.2) exactly converges 
to the function \ipn{r, A)), if when taking sum by m a totality (infinite number) 
of functions \^rn{r}) is considered. However the wave function \tpn{r, A)) from a 
qualitative point of view may be with a high degree of accuracy described by the 
expansion like (3.2) using a finite number N, {N » 1) of states |$TO(0)-Tlie 
further increase of the number of states N is similar to the tendency of the 
Plank constant to zero h ^ Q and leads to the improvement of quantitative fit 
between classical and quantum chaos [4, 5] . 

According to the premises one may say, that the problem of finding eigen- 
functions and eigenvalues of Hamiltonian H{p, f, A) may be solved by a numeric 
diagonalization of the matrix equation (3.4). 

But here a question is born. If the state of the system is described by a finite 
superposition of regular states |<I>m(^)), then the state \ipn{r, A)) is regular too, 
so where should one look for a nonregular behavior of the system? 

By the assumption taken for today [1,28,29] the properties of chaos in quan- 
tum system are revealed in that sense that matrix (3.5) is a random matrix, 
all the elements of which are random values. Based on this, there must exist 
some correspondence between matrix (3.5) and the quantities describing clas- 
sical chaos. Using our problem as a pattern we shall start to examine this 
question. 
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Let us start with quantum mechanical investigation of the Hamihonian (2.2). 
As it has been already mentioned, when the external magnetic field is absentA = 
the Hamiltonian H{p, r, 0) is exactly integrable. When the external magnetic 
field is turned on A > 1 the Hamiltonian becomes nonintcgrablc, and for some 
definite value A = Aq the solutions of the classical equations corresponding to 
the Hamiltonian H{'p, r, Aq) display chaos. By analogy with classical chaos, 
when a set of phase trajectories with a small dispersion of initial conditions 
is investigated, we shall consider a small increment of the parameter 6\ = 
A — Aq. Using the basis of functions in which Hamiltonian H{p, f, Aq) is diagonal, 
for small values of 5X the Hamiltonian H{p, f, A) may be represented as a sum 
of two matrices: 

H = Eo + 6XB, (3.6) 

where Eq is a diagonal matrix, i? is a banded matrix, the elements of which 
are random numbers. By a numerical diagonalization, one can define the eigen- 
functions \tpn{r, A + Aq)), \tpn{r, A)), and eigenvalues -E„(Ao + i5Ao), En{Xo) cor- 
responding to the H{p,r,X) and H{p,r,Xo) respectively. We are interested in 
the interval of energy 6E = E{Xo + SXq) — E{Xo) which is small from a classical 
point of view, and large from a quantum mechanical one (with relation to the 
large number of quantum levels contained in this interval). 

Let us begin to the quantum mechanical analysis of the Hamiltonian (2.1). 
The eigenvectors of the basis, in which Hamiltonian is diagonal, can be found 
easily and are of the form [23] 

|$n,i,m(r, e, if)) = \Rn,iYim{e, ^)) , (3.7) 

where R^i = „-+2(2;+i), V^^^j^(2r)'e-'-/'f (-n + 1 + 1,21 + 2,^) is the nor- 
malized radial part of the wave function and Yim{0, F{—n + 1 + 1,21 + 2, ^) 
are the spherical and hypcr-gcometrical functions respectively [27]; r, 6 and (p 
denotes the spherical coordinates. 

As was mentioned above, using of the symmetric Kepler problem's eigenfunc- 
tions is connected to the requirements of the random matrix theory. Calculation 
of non-diagonal matrix elements in this basis is more complicated then in Lan- 
dau basis [15,23]. But the fact that we are interested in the state I = 0, m = 
makes life and it's still possible to obtain some analytical results. Actually the 
problem is reduced to the calculation of the following matrix elements: 

(*«,o,o(r,6',99)|.x2 +y^|$m,o,o) = ($„,o,o(?-,6',¥')k^sm2^|$„,o,o(r,6',(y5)). (3.8) 

These matrix elements can be calculated analytically. As a result for diagonal 
matrix elements we get [25]: 

2 

($„,o,o(r, e, (p)\r^sin^e\<i>n.o,o) = y^^^^ + 1) (3-9) 

and non-diagonal matrix elements can be reduced to the calculation of the 
following integral: 

(*„,o,o(r, 0, ip)\r^sinH\^mMr, 0, ^)) = 3/2 -^7 '°(^' (3-10) 
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/ k + k I I 

J^'°{a,a) = / e ^ F{a,^,kr)F{a' ,j,k'r)dr, 
Jo 

where 

2/2 

k = —;k = — ■,a = —n + l;a = — m + 1; 7 = 2; 5 = 3. 

n m 

Direct calculation of the integral J^'°{a,a ) is not possible. But we can use 
recurrence formula, which connects integral J^'^{a,a ) and integrals with the 
lower values of superscript J^~^'^{a,a ); J^~^-^{a,a ). After straightforward 
but laborious calculations one can obtain J^'°{a,a ): 

"^7'° = fc2^fc'2 [(|(fc - k) -ka + k'a -k'{S- 1)) Jf -^'°(a, a')+ 

+{S - 1)(7 - 1 + 5 - 1 - 2a') J7 + 2a' (5 - l)J^-^'\a, a' + 1)] 

where 

J0'0(a,a ) = 2'^r(7)(fc + fc -''{k -k)-%k-k)-^ F{a,a 

(3.11) 

and r(7) is the Euler gamma function [27]. 

4 Irreversible evolution of the quantum chaos. 
Formation of the mixed state. 

At the discussion of the nonreversible evolution of the quantum -mechanical 

system, naturally the question emerges. How in a quantum system originates 
non -reversibility? If the quantum system evolves according to the Schrodinger 
equation, how can a pure quantum -mechanical state become a mixed one? The 
question is that, in contrast to the classical chaos, quantum mechanically irregu- 
lar motion cannot be characterized by extreme sensitivity to tiny changes of the 
initial data. Due to the unitarity of the quantum dynamics, the overlap of two 
wave functions remains time-independent |(?9(t)|C(i))P = l(*^(0)IC(0))P, pro- 
vided time-dependence of -dit) and C,{t) is generated by the same Hamiltonian. 
However, an alternative characterization of classical chaos, extreme sensitivity 
to slight changes of the dynamics does carry over into quantum mechanics. 

Let us assume, that the amplitude of the external magnetic field is modulated 
by varying field 

B = Bq + ABocosQt. (4.1) 

Taking into account (4.1) the Schrodinger equation for the wave function of 
diamagnetic hydrogen atom takes the form 
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where 

H{t) = Ho{p, f, Ao) + AXit)V{r, 0), 
AX{t) = AXcosQt. 

The solution of the time-dependent Schrodinger equation formally can be 
written with the help of a time-ordered exponential 

U{t) = {exp[-^J^dtH{t)]}+, 
where the positive time ordering requires 



[A{t)B{t )]+ = 



A(t)B{t) ift>t 
B{t)A{t) iit<t 



In our case H{t + "^jr) = H{t), the evolution operator referring to one period 
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Tq = the so-called Flouqet operator U{Tq) = F [1], is worthy of considera- 
tion since it yields a stroboscopic view of the dynamics 

v;(fcTo) = {Prm- (4-3) 

The Floquet operator, being unitary, has unimodular eigenvalues. Suppose we 
can find eigenvectors |<^^) of the Floquet operator 

F\v^) = e-^^x|^^), (4.4) 
{"PxWe) = he- 

Then, with the eigenvalue problem solved, the stroboscopic dynamics can be 
written out explicitly [1] 

V(fcTo) = J2^-''^H^xmomx)- (4-5) 



As it was mentioned above, our aim is to proof that one of the signs of the 
emergence of quantum chaos is a formation of the mixed state. Being initially 
in a pure quantum mechanical state, described by the wave function tpn, the 
system during the evolution makes an irreversible transition to the mixed state. 

After the formation of a mixed state the quantum-mechanical description 
loses its sense and we have to use a quantum-statistical interpretation. To prove 
the formation of the mixed state, one has to show the zeroing of non-diagonal 
elements of density matrix, the equality of which to zero is a sign of a mixed 
state [30,31]. 

Taking (4.5) into account, we obtain for non diagonal matrix elements of the 
density matrix: 

Pnm = (4.6) 
where Cn = e~*'^^"(0„|V'(O)), (...) means averaging over time. 
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For averaging let us recollect that value 0„ is the eigenvalue of the Floquet 
operator. Owing to the non-integrability of the system, eigenvalues can be 
obtained only by way of numerical diagonalization of the Hamiltonian H (4.2). 
According to the main hypothesis of the random matrix theory [1,11], elements 
of this matrix in the chaotic domain are random values. So, their eigenvalues 
also can be considered as random values. This statement is valid only in the 
chaotic domain [32]. According to this, the phase 

f{n,m) = (t)m- (t>n, (4.7) 
also may be considered as random. Therefore when T = kTo, fc 1, we get 

rr. T/To 

Pnm = Yll (e-*'=^("•'"))(.^„^|V(0))(^(0)|<^„) = 0, (4.8) 

k=l 

where (...) denotes ensemble average of random matrices, that corresponds to a 
small dispersion of the value of magnetic field parameter ^Aq. 

Expression (4.8) may be found by more vigorous mathematical substantia- 
tion. For this let us recall some details from the probability theory [33]. 

a) In general case, under the characteristic function of the random value X, 
mathematical expectation of the following exponent is meant 

F{t) = M{exp{itX)), (4.9) 

where f is a real parameter. 

b) Mathematical expectation itself is defined as a first moment /xi of the 
random variable X 

M{X) = (X) = Ml = ^ XkPk, (4.10) 

k 

where X is the discrete random variable which takes possible values x\,X'2,--- 
with appropriate probabilities Pk = P{X = Xk), (...) means average. 

Taking (4.9), (4.10) into account and considering /(n, m) as a random value, 
we get: 

rj. T/To 

Pnm = YT. nk), (4.11) 
k=l 

where 

F{k) = M(e*'=^("''")), (4.12) 

is the mathematical expectation of the characteristic function of the random 
phase f{n,m), and 

A = {ci,rn\m)mmn). 

Let us assume that random phase f{n,m) has a normal dispersion [33]. Then 
from (4.12) one can obtain 

i^(fc) =e*"'=e-'^, (4.13) 
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where 

o = M(/(n,m)), (4.14) 
is the mathematical expectation of the random phase, 

= M{f{n, m)) - (M(/(n, m)))\ (4.15) 

is the mean square deviation. 

Substituting (4.13) in (4.11) and making transition from summation to the 
integration J2k J dk finally we get 

p„.(r) = ^(i^r/^I-^] +i^r/i[^±ffi^]), (4.16) 

where Co = iAe^ and Erfi{...) is the error function, connected to the ordinary 
error function [33] by the relation 

Erfi{x) = -i^Erfi{ix). (4.17) 

Taking into consideration asymptotical form of the function Re{Erfi{x)) 0, 
X ^ CO and Im{Erfi{x)) ^ 1, x co, in case when distribution of the eigen- 
values of the quantum-mechanical chaotic system obeys the normal distribution, 
we have 

|p_(T)|~l/T~0, 

r>ro. ^ ' ' 

Zeroing of the non-diagonal elements of density matrix is the sign of quantum 
chaos beginnings and formation of mixed state. From that moment the quantum 
dynamics is non-reversible since the information about the wave functions phase 
is lost. Given above reasoning can be checked by numerical experiment, which 
consists in: the direct solving of non-stationary Schrodinger equation (4.2), 
defining eigenvalues (p^ and in estimation of the averaged non-diagonal elements 
of the density matrix according to (4.8). 

The result of numerical calculations is represented on Fig. 6. 



5 Quantum Poincare Recurrences 

According to the results of the previous section (see Fig. 6), for times t > Tc for- 
mation of a mixed state happens due to time evolution of quantum-mechanical 
system. Consequently the concept of a wave function loses its meaning and the 
Poincare recurrences are out of the question. But what happens when t < t^? 
We shall clear up this question in the given section. 

As it was mentioned above, in work [22] the concept of quantum Poincare 
recurrence was introduced. In particular, the authors have studied the quan- 
tum dynamics of one-dimensional rotator under the influence of delta- like per- 
turbation. The definition of the "distance" between the wave- functions was 
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Figure 6: The graph of the dependence of the second term pmn of non-diagonal matrix 
element (4.6) on time t = kTo. The graph is obtained using (4.8) and by way of nu- 
merical diagonalization of the Hamiltonian (4.2) for different values H{ti) = H{X{ti)), 
ti G [0,T]. The graph is plotted after taking an ensemble average of Hamiltonians 
H{Xo), corresponding to the dispersion 5Xq = 0.1 of the value of magnetic field Ao. 
As it is easily seen from the graph that when t = kTo > Tc, To — ^ — O.Itc, zeroing 
of Pnm happens. In the numeric calculations we have used definite initial conditions 
tpn{0), V'm(O), obtained via numerical diagonalization of the Hamiltonian (2.1) for 
Ao = 1. The dimension of the diagonalized matrix was 500 x 500, Xj C = Ij 500 

introduced 

(f{t) = {^{r,t) - 4'(r,0)|*(r,i) - *(f,0)), (5.1) 
where |^'(r, t)) is the solution of nonstationary Schrodinger equation 

a* 

ih— = F(t)*. (5.2) 

The authors of the work, cited above, look for the solution of the equation (5.2) 
in the form of expansion of the function \'^{f,t)) in an arbitrary orthogonal 
basis |$to) 

\'if{f,t))n = ^Cnm{mmm}, (5.3) 
m 

where Cnm{t) are the expansion coefficients, which determine the time depen- 
dence of the function |^'(r, It is obvious, that during the cvohition in time 
of the wave-function |5'(f, t)), the quantity d{t) will change in time. The au- 
thors offer to consider the quantum Poincare recurrences as such values of time 
t, that satisfy the condition \d{t)\ < e, where e is a small constant value (e < 1). 

Our purpose is the use of this method for the study of quantum Poincare 
recurrences in the system of diamagnetic hydrogen atom. For that, like in the 
previous section, we shall assume, that the amplitude of the external magnetic 
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Figure 7: The histogram of distributions of the times of quantum Poincarc recur- 
rences, plotted in compliance with the definition (5.1). The quantity N{t) determines 
the dependence of the number of quantum Poincare recurrences on the distribution of 
the periods of recurrences. Ao = 1, AA = 0.1, H = §2:, To ~ O.Itc « 0.5 • 10"^^s. 



field is modulated by a varying field. But in contrast to the previous section, 
we shall examine the solution of the equation in interval t < Tc, when the 
system is still in a pure quantum-mechanical state and chaotization of phase 
and formation of a mixed state can be neglected. 

Before starting the analysis of the equation (4.2), let us note some typical 
differences of our problem from the problem studied in [22] . In problem studied 
in [22], the dimension of the Hilbcrt space is finite. That is why the wave- 
function |\l/(r, t))„, according to (5.3), is exactly defined by the finite set of 
the orthogonal functions The authors make an assumption, that the 

expansion coefficients Cnm{t) are the random quantities and satisiy the Gaussian 
distribution. 

In our case the dimension of the Hilbert space is infinite and the approxima- 
tion of i))„ by the finite set of functions of the orthogonal basis |<l>m(r)) 
is rough. Our purpose is a proof that in spite of it the quantum Poincare 
recurrences may happen in the system. 

For the solution of the nonstationary Schrodinger equation (4.2) we shall 
consider i as a discrete quantity (consequently the values of the parameter 5\{t) 
become discrete too). We shall not assume the coefficients C„m(A,t) to be 
a priori random, but we shall determine their values by solving the matrix 
equation (3.5) and taking into account the dispersion of the parameter 5\o. 
Thus, using the formulas (3.2), (3.5) and (5.1) one can define the distribution 
of the times of Poincarc recurrences. The results of numeric calculations are 
shown in Fig. 7 and Fig. 8. 

According to Fig. 8, the distribution of the times of Poincare recurrences is 
of the Lorenz form. This fact can be easily explained, taking into account the 
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Figure 8: Exponential approximation of tlie numeric data represented in Fig.7, made 
in compliance with the formula (2.4). The graph is plotted in a logarithmic scale. 

Therefore a straight line on the given graph corresponds to the exponential Lorontz 
form. The result allows to determine numerically the mean time of quantum Poincare 
recurrences and to compare it with the mean time of classical Poincare recurrences 
Tree ~ T ~ 1 //lo (see formula (2.4)). As it comes from the graph, ~ ~ 

0.5 • 10~^'*s, where 6 is the angular coefficient of the approximating straight-line. The 
achieved estimate for the mean time of quantum recurrences is in a good agreement 
with the mean value of the classical Poincare recurrences Tree ~ "^c ~ ^ = 0.5 • 10~^*s 
(see Fig. 3). 

existence of a strong classical chaos (i.e. the absence of the stability islands in 
the phase space). 

The Lorentz form testifies that \i t > = Tc, quantum Poinacre recur- 
rences are improbable. This coincides with the results of the previous section, 
according to which the system is found in a mixed state in this time interval. 

Based on the results obtained above, one may assert that quantum Poincare 
recurrences also take place in quantum systems with infinite "phase space" (i.e. 
with infinite-dimensional Hilbert space). 
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